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Abstract 

Comparative studies of magnetoelectric susceptibility (a), magnetization (M), and magnetostric- 
tion (u) in TbMn205 reveal that the increment of M owing to the field-induced Tb 3+ spin align- 
ment coins a field-asymmetric line shape in the a(H) curve, being conspicuous in a low temperature 
incommensurate phase but persistently subsisting in the entire ferroelectric phase. Correlations 
among electric polarization, u, and M 2 variation represent linear relationships, unambiguously 
showing the significant role of Tb magnetoelastic effects on the low field magnetoelectric phenom- 
ena of TbM^Os. An effective free energy capturing the observed experimental features is also 
suggested. 
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Nontrivial cross- coupling between electric and magnetic dipoles realized in multiferroics 
has been a subject of extensive researches in recent years, which are targeted to understand 
the mechanism of magnetoelectric (ME) coupling as well as to find novel device applications. 
One of the key compounds that has triggered such research activity is TbMn 2 05, in which 
a continuous actuation of electric polarization (P) is realized within low magnetic field (H) 
below 2 T. Numerous studies on this compound and related -RMn 2 Os (R = Y, Dy, Ho, 
Er, and Bi) have shown that spatially modulating, noncollinear magnetic order due to spin 
frustration is responsible for inducing ferroelectric order in these materials. More specifically, 
a main mechanism for having nontrivial P in i?Mn 2 Os (R = Tb, Y, Dy, and Bi) has been 
attributed to exchange striction among frustrated Mn spin networks,- - - while P contribution 
from spiral spin order has also been known to be important in i?Mn 2 Os (R = Ho, Er, and 
Tm). 6,7 Thus, a main mechanism for developing P in i?Mn 2 Os can be arguably dependent 
on a specific material while it is obviously associated with the Mn spin order.- - - 

Only a limited number of works have discussed the possible effects of rare earth ions on 
the temperature- and if -dependence of P on i?Mn 2 Os.-~— As a result, a proper role of rare 
earth ions on the ME phenomena of i?Mn 2 Os is far from complete understanding and thus 
worthy of investigation. One particularly important question is how one can understand the 
H- induced actuation of P that is uniquely realized in TbMn 2 5 . A detailed understanding 
of this intriguing question is likely to provide not only an answer for the long-standing 
puzzle that has triggered the multiferroic research but also useful information regarding the 
application of multiferroics. 

In this communication, on the basis of systematic studies of magnetostriction (u), mag- 
netization (M), and ME susceptibility (a), we uncover that M change due to Tb spin 
alignment with H determines the evolution of both u and P predominantly, thereby devel- 
oping linear relationships among M 2 , u, and P in the entire ferroelectric phase. An effective 
free energy analysis based on the magnetoelastic coupling of Tb can successfully describe 
the experimentally found correlation among those physical quantities. 

Single crystals of TbMn 2 Os were grown with a PbO:PbF 2 flux.— To investigate detailed 
if-and temperature-dependent P (/ /b) and lateral length I (//a) change, we developed a 
sensitive ME susceptometer and a high precision dilatometer, both of which operate in a 
PPMS™. In this study, we have focused on a 21 = 5Pb/SH a and longitudinal magnetostric- 
tion u a = {I (Ha) — 1(0))/ 1(0) along the a-axis. For the former, we used solenoid coils to 
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FIG. 1. (color online) (a) Temperature dependence of electric polarization (P) at selected magnetic 
fields. Tq and T\q refer to the commensurate and incommensurate antiferromagnetic transition 
temperatures, (b) Tb spin configuration of TbM^Os at 20 K reproduced from Ref. 4. (c) 
Magnetization (M) along the a-axis at 3, 20, and 40 K. 

apply small ac H (/ /a) of ~4 Oe and a high impedance charge amplifier to sensitively de- 
tect an ac modulated charge, proportional to 5Pb, by using a lock-in technique.— Dielectric 
constant (e) and M were also investigated with a capacitance bridge and a vibrating sample 
magnetometer, respectively. 

Upon cooling, TbMr^Os passes through three main magnetic and electric transitions: 
an incommensurate (ICM) magnetic ordering at Tn ~ 43 K, and a nearly concomitant 
ferroelectric and commensurate (CM) magnetic ordering at Tq « 37 K, and a reentrant 
low temperature incommensurate (LT-ICM) magnetic ordering at Tic ~ 25 K with a sharp 
decrease in P& (See, Fig. l(a)).->^ Those transitions are also accompanied by structural 
anomalies.— The ferroelectricity and the structural instability is postulated to stem from 
atomic displacements of Mn 3+ ions located at the centers of bipyramids.- 1 ^ Although the 
antiparallel alignment of Tb spin moments, as shown in Fig. 1(b), has been extracted 
from the neutron scattering refinement below ~20 K,-^ there is no clear evidence for a 
thermal transition of Tb spin ordering below Tb, in contrast to the case of Dy spins in 
DyMr^Os.- Moreover, the three thermal transition of TbM^Os are quite similar to those 
of an isostructural YMn 2 Os without any rare earth ion.- Owing to these facts, the effect of 
Tb 3+ ions on the physical properties of TbMn 2 Os appears small. 



However, there exist a couple of experimental features that warrant explanations based 
on the Tb spin effects on TbMn 2 C>5. First, in contrast to YMn 2 05, in which the negative P& 
hardly changes up to 9 T in the LT-ICM phase, 16 P& of TbMn 2 05 increases with decreasing 
temperature below ~15 K at HoH = T, and this low temperature positive P& is drastically 
suppressed to become negative at fi H = 2 T and even more at 9 T.— Second, the temper- 
ature range for the Pb increase is consistent with that of the Tb moment increase observed 
by neutron scattering, thereby indicating a nontrivial coupling between Pb and Tb spins.- 
Third, as shown in the isothermal M vs. H curves in Fig. 1(c), the spins of Tb 3+ ions align 
within fi H ~ 2.5 T at 3 K and ~ 8 T even at 20 K to nearly reach a predicted saturated 
moment (M s ) of 9/ie/f.u. (4/ 8 , 7 P 6 ) and thus, the Tb spin alignment is a dominant source 
of MM 

The large M due to the Tb spin alignment results in a significant change in length, 
under H, i.e., magnetostriction in TbMn 2 Os. Figure 2(a) shows that the u a is positive 
and increases in proportion to M 2 . The u a value of +6 x at 2 T is indeed similar to 
the longitudinal magnetostriction observed in compounds with the Tb 3+ ions; for example, 
longitudinal striction is +2 x 10~ 5 in TbA10 3 at 4 T and +5 x 1(T 5 in Tb 3 Ga 5 0i 2 at 2.2 
T.— ^ According to these two features in u a , it is most likely that the magnetostriction of 
TbMn 2 05 is mainly attributed to Tb 3+ ions involving both single as well as two ion effects 
as in the case of TbAlOs.— 

M, u a , and P variation under H is closely linked to the Tb spin moment. The a 2 i(P) 
curves in Fig. 2(b) directly show an evidence for such nontrivial effects of Tb spin moment 
on ME phenomena. a 2 i(P) at 3 K displays a sharp dip and peak structure around ±0.6 T. 
Upon being integrated with H as Pb{H) = P>(0) + J^ 1 a 2 i<iP, Pb{H) at 3 K steeply decreases 
within | if | < 2 T (Fig. 2(c)), which is consistent with the reported data from pyroelectric 
current measurements.— 1 ^ The decreasing Pb{H) turns out to be proportional to M 2 in a low 
H region as is the increase in u a , thereby establishing an unambiguous and close correlation 
between the decrease in Pb and increase in u a at 3 K. This correlation is further corroborated 
by the close similarity in the characteristic asymmetric line shape observed in both -du a /dH 
and a 2 i(if) curves. Although the absolute value of u a is too small to directly account for 
the absolute change of P&, this correlation reflects that Tb-0 distribution can be changed 
by a local strain of Tb 3+ ions~i or exchange interaction between Mn and Tb ions further 
modulate spin ordering patterns of Mn 3+ ions^&i2 to amplify the concomitant Pb decrease 
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FIG. 2. (color online) -ff-dependence of (a) longitudinal magnetostriction u a , (b) magnetoelectric 
susceptibility 021 = 5Pb/5H a , (c) Pb determined from the integration of 021 with H, and (d) 
-dua/dH at 3, 15, and 26 K. The dashed lines in (a) and (c) represent scaled M 2 curves at 3 
and 26 K to fit into the low field. A dashed line in (d) represents a guide to eye to illustrate the 
asymmetric line shape of -du a /dH at 26 K. (e) Enlarged at2i(H) at positive H region at various 
temperatures. Dashed line is a guide to eye. 

under H. All these observations consistently support that the ME phenomena of TbMn 2 05 
at 3 K are coupled with magnetostriction mainly due to the Tb 3+ ion in a nontrivial way. 

It is further noteworthy in Fig. 2 that the magnetostriction effects of Tb 3+ ion seen in 
-du a / dH and a 2 i {H) curves are well maintained up to high temperatures. Except the large 
peaks in the -du a /dH due to the Mn spin transition from the LT-ICM to CM states at 26 
K, the asymmetric line shape of the -du a /dH curves is clearly observable at 15 and 26 K 
(Fig. 2(d)), signaling a significant magnetostriction effect in the entire ferroelectric phases. 
The ot2i{H) curves at 15 and 26 K (Fig. 2(b)) also show the characteristic asymmetric line 
shape, except jumps at 26 K that are coming from the same Mn spin transition. Similar 
to the relationship between and M 2 , the asymmetric line shape of -du a /dH results in 
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FIG. 3. (color online) Temperature (T) vs. magnetic field (H) phase diagram of TbM^Os. 
Asterisks and squares indicate the phase boundaries determined by the a and e measurements, 
respectively. Solid and open symbols represent the data measured during the H- or T-increasing 
and decreasing runs, respectively. Solid circles and diamonds refer to the points of ot2i{H) minima 
and e(H)H maxima, respectively. 

the characteristic increase in u a proportional to M 2 , thereby demonstrating the nontrivial 
coupling between P& and u a at temperatures below Tq. 

To estimate the phase region affected by Tb 3+ magnetostriction, we trace the character- 
istic minimum positions, a2i,min(P), seen in the asymmetric line shape of a2i{H) at H > 
(solid circles in Fig. 2(e)) and plotted in the phase diagram of Fig. 3. The a^minO^) ex ^ a ^ 
all temperature regions below Tq. The phase boundaries for the LT-ICM to CM transitions 
of Mn spins are also determined from the hysteretic jumps in the at2i(H) (asterisks in Fig. 
2(e)), ct2i(T) curves (not shown here), and in our previously published e(T) datar^-. In the 
CM phase region, the trace of the a2i,mm(H) (solid circles) is significantly shifted to higher 
fields, thereby indicating that the complete alignment of Tb spins becomes easier at low 
temperatures due to the increment of thermal entropy in the LT-ICM phase. e(H) showed 
a maximum, e max (H), of which trace was determined from the results in Ref. 17 (diamonds 
in Fig. 3). Although e max (if) is shown at somewhat larger H, it shows a similar curvature 
change as the trace of «2i,mm(-^), thereby indicating that Tb magnetostriction also affects 
the magnetodielectric effect. 

We further uncover that the isothermal variation of P&, u a , and M 2 in a low field region 
roughly follows a simple relationship, i.e., oc u a and M 2 . Figure 4(a)- (c) shows a com- 
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FIG. 4. (color online) Plots of (a) pb vs. u a , (b) pb vs. m 2 , and (c) u a vs. m 2 at 3, 5, 15, and 26 K. 
Dashed lines are linear fit lines for low field region, (d) Comparison of the temperature-dependence 
of —dpbjdu a (left, circles) and —dpb/dm 2 (right, solid squares). 

parison of three unitless quantities pb, u a , and m 2 . Here, for the convenience of description, 
we define pb = Pb/Pm&x and m = M/M s with P max = 42.3 nC/cm 2 (P at 26 K and T) and 
M s = 9yUB/f-u. As shown in the figures, the variation of pb, u a , and m 2 is roughly linear to 
each other, except a jump due to the Mn spin transition. We further note in Fig. 4 that 
there exists characteristic temperature-dependence in their linear relationship; the linear 
slopes of pb vs. u a and pb vs. m 2 curves show strong temperature dependence while those 
of u a vs. m 2 curve is nearly temperature-independent. The overlapping of —dpb/du a and 
—dpb/dm 2 curves at all temperatures below Tq with a single constant multiplication (Fig. 
4(d)) confirms that the slope changes in both p b vs. u a and pb vs. m 2 curves follow almost 
the same temperature-dependence. 

To understand the intriguing coupling among P{= Pb), M, and u{= u a ), we consider a 
free energy that effectively considers the magnetoelastic effect of Tb 3+ ions. 

F H (P, M, u) = (P ~ Ph=o)2 + ^- + \c T u 2 - -P 2 M 2 - f -P 2 u - gM 2 u - ME (1) 

2Xe,T 2Xm,T 2 2 2 

The first three terms describe temperature-dependence of the order parameters P, M, 
and u. To describe the variation of quantities under low if-regime, we assume that Ph=o, 
Xe,T, Xm,T, and Ct have predetermined temperature-dependence, consistent with the exper- 
imental data. The temperature-dependent evolution of Ph=o and Xe,r has been well studied 
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and appears to be mainly determined by the Mn spin interactions,-^^ while Xm can be 
determined by Mn-Mn, Tb-Tb, and Tb-Mn interactions.- However, to our knowledge, elastic 
stiffness constant (CV) has not been known yet. The next three terms describe couplings 
among the order parameters. These are invariant with both inversion and time reversal sym- 
metry operations. Here, A, /, and g correspond to the temperature-independent coupling 
constants that are specific to the material. 

By minimizing Eq. (1) with M, the usual form of M — Xm,TH is obtained under a 
condition of l/xm,T 3> (2gu + \P 2 ). With the replacement of M with Xm,rH, F becomes a 
function of P and u. The simultaneous minimization of F with respect to P and u further 
provides two linear equations with variables P and u. By obtaining the functional form 
of dP/dH 2 and du/dH 2 from the two linear equations, one can represent du/dM 2 and 
dP/dM 2 , as shown below. 

du - fu ~ Xx 2 m , T H 2 )g + f\P 2 du 

(2) 



xl, T dH 2 (^-fu- Xx 2 m , T H 2 )C T - f 2 P 2 dM* 



dP fgP + XC T P dP (3) 



X 2 m , T dH 2 (^-fu- \x 2 m , T H 2 )CT - f 2 P 2 dM 2 
According to the results shown in Figs. 4 (b) and (c), du/dM 2 is temperature- independent, 
while dP/dM 2 is not. To satisfy both of these constraints, the second term in the de- 
nominator might be dominant over the first, which then results in du/dM 2 = —A//, 
dP/dM 2 = -(fg + XC T )/f 2 P, and finally, dP/du = (fg + XC T )/fXP. In this approx- 
imation, Eq. (1) can successfully explain the temperature-independence of du/dM 2 as well 
as the same temperature-dependent variation of dP/dM 2 and dP/du since the Ct and P 
commonly determine the temperature dependence of the last variables. Therefore, our ap- 
proach of an effective free energy, which is based on the magnetoelastic effect of Tb spins, 
can provide explanations on the intriguing coupling relationship among P, u, and M. 

All above experimental results and considerations based on the free energy constitute 
compelling evidences for the existence of significant magnetoelastic effect due to Tb 3+ ion 
in TbMn 2 05 so as to modulate the macroscopic physical quantities of M, P, and u, which 
correspond to the spin, charge, and lattice degrees of freedom in TbMn 2 5 , respectively. 
The ME phenomena of TbM^Os such as the actuation of P and the variation of P propor- 
tional to M 2 can only be explained by the magnetoelastic effects of Tb spin moment. Our 
results provide several implications for the physics of TbMn 2 Os and related multiferroics. 
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First, temperature-dependent elastic constant Ct can be directly linked to the temperature- 
dependent variation of dP/dM 2 and dP/du via above Eqs. (l)-(3). This scenario can be 
checked from a direct measurement of Ct- Second, a microscopic mechanism of how the 
Tb spin alignment with H can amplify the P variation is a subject worthy of further explo- 
ration. The issue can be linked to either a local strain field of Tb magnetostriction or direct 
exchange coupling between Mn and Tb spins. Third, for a proper description of multiferroic 
phenomena as well as its application, consideration of the magnetoelastic effects of magnetic 
ions can be generally important in many other multiferroic materials as well.— 

In conclusion, we have provided clear experimental evidences and theoretical indications 
that magnetostriction due to the Tb spin alignment crucially affects the ME phenomena of 
TbMn 2 05 in the entire ferroelectric phase. Our results imply that a proper control of the 
strain or magnetic moment of rare earth ions can be useful in the application of existing 
multiferroics in a low field phase. 
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